
STOP 



Early Journal Content on JSTOR, Free to Anyone in the World 

This article is one of nearly 500,000 scholarly works digitized and made freely available to everyone in 
the world by JSTOR. 

Known as the Early Journal Content, this set of works include research articles, news, letters, and other 
writings published in more than 200 of the oldest leading academic journals. The works date from the 
mid-seventeenth to the early twentieth centuries. 

We encourage people to read and share the Early Journal Content openly and to tell others that this 
resource exists. People may post this content online or redistribute in any way for non-commercial 
purposes. 

Read more about Early Journal Content at http://about.jstor.org/participate-jstor/individuals/early- 
journal-content . 



JSTOR is a digital library of academic journals, books, and primary source objects. JSTOR helps people 
discover, use, and build upon a wide range of content through a powerful research and teaching 
platform, and preserves this content for future generations. JSTOR is part of ITHAKA, a not-for-profit 
organization that also includes Ithaka S+R and Portico. For more information about JSTOR, please 
contact support@jstor.org. 



Investigations in the Lunar Theory. 

By Ernest W. Brown, 
Professor of Applied Mathematics at Haverford College, Pa. 



This paper is an outline of a plan for the development of the expressions 
which represent the coordinates of the Moon, together with certain theorems 
connected with the infinite determinants which determine the motions of Perigee 
and Node, and, in addition, some results concerning the constant part of the 
expression which gives the Parallax of the Moon. 

It has been pointed out that the algebraical expressions of the coordinates 
show slow convergence chiefly when the coefficients of the periodic terms are 
arranged in powers of the ratio of the mean motions of the Sun and Moon, and 
that when we use the numerical value of this ratio from the start, keeping the 
other constants arbitrary, slow convergence in the series arranged in powers 
and products of the other constants, is unusual; also, the observed value of 
the ratio of the mean motions is determined with an accuracy far surpassing that 
with which the other constants are known. These reasons have seemed to point 
towards a semi-algebraical development, in which the ratio of the mean motions 
is given its numerical value while the other constants are left arbitrary. 

The developments given below are built up on this basis. The orbit of 
Dr. Hill, which depends on on only, is, therefore, a numerical one, and is used as 
a first approximation or " intermediate " orbit. 

But there are several methods by which we might proceed with the develop- 
ment from this point. It is possible to use polar or rectangular coordinates, and 
the equations of motion in each case can be given several different forms. Of 
the forms which the equations may take when rectangular coordinates are used, 
there are two which seem specially adapted for continued approximation. Of 
these one set are of the second degree and second order when we neglect the 
Parallax of the Sun and the Latitude of the Moon. "When these quantities are 
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included, two equations still keep a form in which the principal terms are of the 
same kind, while the third equation is a sufficiently simple one. These are the 
equations (8) below. Some of the elliptic and parallactic terms have been 
determined from these equations. 

Another method is to use the equations of motion in their oi'iginal form. 
This requires the expansion of such expressions as (x + Sx)/(r + Srf, where 
x , r represent the values of a;, r in the intermediate orbit. These being trigo- 
nometrical or exponential series, entail a certain amount of labor in the develop- 
ments. How this is effected and in what way we gain or lose in comparison to 
the former method is shown with some detail in Part I. It is on these lines that 
the developments there given are carried out. 

The chief difficulties which arise in any method of treating the Lunar Theory 
are those connected with the determinations of the motions of the Perigee and 
Node. When Hill and Adams conceived the idea of the infinite determinant, and 
succeeded in solving one, there appeared to be an opening for the determination 
of these motions with an accuracy which should leave nothing to be desired. 
The equation for the latitude in rectangular coordinates is in a form which gives 
the principal part of the motion of the Node — that depending on m only — by 
means of an infinite determinant directly, and thence the coefficients of the 
periodic terms are easily obtained. But in the case of the Perigee this is not so. 
Several transformations are necessary before a suitable determinant is evolved, 
and when these have been made and the coefficients of the transformed solution 
obtained, the return process to find the coefficients of the periodic terms in 
rectangular coordinates is a laborious one. Once the motion of the Perigee has 
been obtained, however, we can use the equations in any form we choose. 

So far the question, then, of finding the principal parts of the motion of 
Perigee and Node is resolved, and the coefficients dependent on these can be 
obtained. But we are again met by the same difficulties when we try to find 
further approximations to these motions. We again obtain an infinite set of 
linear equations. But they now contain constant terms, and though an infinite 
determinant can be obtained for the new part of the Perigee or Node to be found, 
it is in a form which does not admit of easy solution. The equations can be 
solved by continued approximation, but this is very laborious. The difficulty 
has been turned by a slight artifice which eliminates all the unknown quantities 
from the set of equations but that required, and the resulting equation is one 
which does not present any unusual difficulties in calculation. 
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The equations thus obtained not only showed methods by which the motions 
of Perigee and Node can be obtained at any stage by simple processes, but 
indicate the way in which Adams' theorems connecting these motions with the 
constant part of the expression of the Parallax arise. The remarkable proof 
Adams gives of his theorems is well adapted to the end in view, but it shows no 
indication of any possible extension to the higher terms. The methods given 
below, though cumbrous, have perhaps the advantage of shedding a clearer light 
on the theorems and of giving them a form which admits of further development. 
The investigation of these occupies Part II. Other theorems are also proved 
which may be of value when a verification of results is required. 

Part III contains one or two deductions from the infinite determinant con- 
cerning the orders of the coefficients, and especially deals with the cases in 
which short-period terms may have large coefficients. 

For convenience, the papers referred to in the following pages — 

On the Lunar Theory : 

I. G. W. Hill. — Researches in the Lunar Theory. American Journal of 
Mathematics, Vol. I, pp. 5-26, 129-147, 245-260. 

II. Gr. W. Hill. — On the Part of the Motion of the Lunar Perigee, etc. Acta 
Mathematica, Vol. VIII, pp. 1-35. 

III. J. C. Adams. — On the Motion of the Moon's Node in the case, etc. 
Monthly Notices R. A. S., Vol. XXXVIII, pp. 43-49. 

IV. J. 0. Adams. — Note on a Remarkable Property of the Analytical 
Expression, etc. Monthly Notices R. A. S., Vol. XXXVIII, pp. 460-472. 

V. C. Delaunay. — Theorie de la Lune. M6moires de l'Academie des Sciences, 
Vols. XXVIII, XXIX. 

VI. C. Delaunay. — Note sur les mouvements du perigee et du noeud de la 
Lune. Comptes Rendus, Vol. LXXIV, pp. 17-21. 

VII. E. W. Brown. — The Parallactic Inequalities in the Lunar Theory. 
American Jour. Math., Vol. XIV, pp. 141-160. 

VIII. E. W. Brown. — The Elliptic Inequalities in the Lunar Theory. 
Amer. Jour. Math., Vol. XV, pp. 243-263, 321-336. 

And on the Infinite Determinant : 

IX. H. Poincare. — Sur les determinants d'ordre infini. Bulletin de la 
Societe math, de France, Vol. XIV, pp. 77-90. 
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X. Helge von Kock. — Sur les determinants infinis, etc. Acta Math., Vol. 
XVI, pp. 217-295. 

— are gathered together and the reference is made by means of the number in 
Roman numerals standing before each. 

The portion of I dealing with the intermediate orbit is also found in 
Tisserand's Mecanique Celeste, Vol. Ill, Chapter XIV. The Memoir II is given 
in Chapter XV, and III and IV in Chapter XVI of the same volume. 

PAET I. 
Development op the Theory. 

1. The Differential Equations. 
The notations to be adopted below are as follows : 

n , n', the mean angular motions of the Sun and Moon about the Earth ; 

2a', r 1 , e', the major axis, radius vector and eccentricity of the Sun's orbit, sup- 
posed elliptic ; 

x, y, 2, the coordinates of the Moon referred to rectangular axes, of which those 
of a;, y are moving in the plane of the Sun's orbit with angular velocity 
n 1 , and that of z is perpendicular to this plane. The positive direction 
of the a;-axis is directed to the mean place of the Sun ; 

v , the solar equation of the centre ; 

S = x cos v + y sin v ; 

u = x + y*S — 1 , s = » — V*/ — 1 > r s = «s = jc 2 + ?/ 2 ; 

n'' a'*, the mass of the Sun ; 

(i, the sum of the masses of the Earth and Moon in the same units ; 

v = n — n', m = n/v, x = (i/v i ; 

y—. ( ,v{t-h)V=I n—yd ■__ 1 d 

C,-e , ^_^_^__ df . 

With the limitations here adopted — the same as those considered by Delaunay 
(V) and Hill (I) — the potential function is 



(i , n' a' n' a' „ 

11 - Vr 2 + a 2 -v/|V — 2Sr , + r*+z*] /* ' 



41 
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and the equations of motion are 



x 



o / • ,» 3ft 

— 2n'y — n x = -~~ , 
ox 



y + 2n'x — n'*y = -^— , 

3ft 

ft' = ft + *n" (<c 2 + y % ) = ft + *n'V, 
the equations may be written 

X - 2n V=-dx-' 



If 



y + 2n'x 



_3ft' 

~ dy ' 

_3ft' 

3z • 



Transforming to the independent variables u, s, z and the dependent £, 
ttions become 



(1) 



the 



equations 



with the expression 

ft'=-J 



D 2 m + 2mDu = 9- . -~— , 

?> a 3s 

D 2 s - 2mDs - — -4- • 23- , 
v 3 3w 

JL ^ 

r 2 ' 3z ' 



Z> s z 



Vws + s 2 ^ Vr' 2 — 2£r' + ws + z 2 / ^ 
This becomes by expansion 

4- ft' = , ^ , 8 + m 8 ^ + m 2 i [3# 2 — «s - <] 
^ 2 vVus + z* r L J 

+ J?£. ^1 [5£ 8 — 3tf(w + a 2 )] +. . . . 

VMS ■+■ 3^ 

where 

ft, = 3m 2 [-^ £ 2 - i (w + s) 2 ] — w 2 (ms + s 2 )(-^- - l) 

<m 3 /7.' 4 „ 



+ -^ • ?r [.5S 3 — SS{us + * 2 )] + . . . . (2) 
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It is necessary, before going further, to see in what way Q. x involves u, s. 
We have 

S* = (as cos v -+- y sin vf 



= i (u + s) 2 — | (w 2 + s 2 ) sin 2 v + 



4V— 1 



(w 2 — s 8 ) sin 2u, 



and if 






the first two terms of Hx become 
_1+P 



3m 2 [£. ( M + s) 2 — i±£( M 2 + . s 2 ) sin 2 « + 1+ £ (m 2 — s 2 ) sin 2w"l— pm 2 (tts+z 2 ). 

From the known properties of elliptic motion p , sin u are of the order e 1 at least. 

Hence fix is of the order e' or 1/a' at least. When we neglect If a', Hi takes 

the form 

Ox = i (^w a + 25m + <7* 2 — 2proV) , (3) 

where A, B, G depend on e 1 and the angle n't + e' only. When 1/a' is not 

neglected we can put 

&i = Ojj + % -f- (3') 

where o g is a homogeneous function of u,.s, z, of degree q free from fractions. 
Also since the time enters explicitly into &! only through the coordinates of the 
Sun, when e 1 = 0, £l t does not contain the time explicitly. 
In the general case the equations may now be written 

_ 3^i 



(Z> 2 + ImD + f w 2 ) u + |m 3 s- 
(Z> 2 — 2mD +fm> + fm 8 M' 
(D 2 — m 2 )z 



(ms + z 2 ) 1 3s ' 

xs d£i x 

(ws + z^f ~~ 3w ' 

xz 1 8nx 

' (us+zj ~ * dz ' 



(4) 



which are the equations to be used in developing this theory. 

Taking the general value of £l 1} multiply the first of these equations by s, 
the second by u and subtract, we obtain 



D (uDs — sDu — 2mus) +fm J (« 2 -s J ) = «^-«^ 



ds 



du 



(5) 
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With the same multipliers for the first two equations and the multiplier 2z for the 
last, we have by addition 

sD % u + uD*s + 2zDh -f- 2m (sDu — uDs) + |»i 8 (m + sf — 2mV 

2« / arii . an, 

" dz J 



/ dl\ . 312, , 3I2A 
/ , 2v = — ( s -g- 1 + M -Q- + z -rr" ' 



J- (6) 



=—._>>* b ^( 3 ')- 



5 = 8 



Again, multiplying the three equations by Ds, Du, 2Dz respectively and adding, 
we obtain 



D [Du. Ds + (Dzf + f m 2 (u + s) 2 — wV + , 2 * ,., 1 

L ' (us 4- z'YJ 



(us -\- z 2 ) 



an. 



l^+f^ + f^] 



Now in the general case £1^ is expressed explicitly in terms of u, s, z, t, hence 

D£l 1 =^Du + ?2l.Zfc + ?pl Dz + ^1^, 

3m 3s az dt 

and therefore the right-hand side of the previous equation is 



3A 

dt 



Dt — DHj 



3^ 



x^ 



rV — 1 * dt 

= z) , n 1 — dd.! 

= D[D- 1 (z>'n 1 )-n 1 ], 

where D'Hx denotes the operation D performed on £l x only with reference to t 
(or £) so far as it occurs explicitly in £l 1) and D~ x denotes the operation inverse 
to D, i. e. integration with respect to £ followed by a division by £. 

If, then, with these substitutions we integrate, the equation becomes 

Da.Ds + {Dzf + f m 8 (« + sf - mV + ^^^ =0 — 0, + D'^D'O,). (7) 

Adding this to (6) and writing down again equations (5) and the last of equations 
(4), the three transformed equations of motion are 

D 2 (us + z*) — Du.Ds— (Dzf— 2m(uDs— sDu) + im*(u +sf— Zm*z % 



= <?' — £ (? + X K + # _1 0^0 



4 = 2 



3n i 



31^ 



Z> (wZte — sZ)m — 2mws) + f m 2 (« 2 — s 2 ) = s -^ 'w -^— > 

27V 



wi 2 s 



xa 



(ws -4- s 2 )* 



3r ~37' 



(8) 
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There are some special cases to be noticed : 

(i) d = 0, the last term of the first equation disappears and g> 2 = ; 
(ii) — = 0, then o g = 0, except for the value q = 2, £i x = w 2 ; 
(iii) 2 = , the third equation disappears ; 

(iv) e' = 0, —j =0, all the terms on the right-hand sides disappear except 

the constant O ; 

(v) e 1 = 0, l/a' = 0, z = 0, the equations reduce to those studied by 
Dr. Hill (I, II) and they give those inequalities which depend one,m only. 

We are going chiefly to study the equations of motion under the form (4), 
and it will be advisable to write down the special forms to which they 
reduce for zero values of certain of the constants. In the majority of the cases 
to be treated 1/a' is neglected, so that £l x has the value (3) and the equations 
become 

(D a + <ZmD + f m % ) u + f mh + Bu + Gs - , + ^ = 0, 

(It is not necessary to write down the s-equation, since u is a complex quantity 
and the ^-equation includes both the a-equation and the ^-equation.) The par- 
ticular cases of (9) are 

%u 
2 = 0, (Z> 2 + 2mD + |m 2 ) u + f m»* + Bu + Gs — ^ = 0, 



(9) 



9(i) 



g = 0, </ = 0, 



<=' — 



o, 



_B = 0, (7 = , p = in equation 9 (i), 
B = 0, (7=0, p = in equations (9). 



9(H) 
9 (iii) 



We shall require the equations (9) expressed in terms of x, y, z. They are 
x — 2riy — Zn ,% x + A'x + B'y= — p ^y f , 

+ jy«+ ^ = - (ya ^ g2)f , 



y + 2n'jc 

2 + w' s 2 + n' a pz 



(r a + z a ) 



,3\i> 



(10) 
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where A!, B', 0, p are quantities depending on e 1 and arising from the part a 2 of 
the disturbing function. 

2. Method of Development. 

We take first the equation 9 (ii). It is known (Hill, I) that in this case the 
equations admit of a solution which can be expressed in the form 

w = w = a 2a i £ 2i+1 



s 



= u — a Xa^ u + l 1 . , . 

= ,„=a 2a_ i _ 1 ^ +1 J -....+•» (11) 



where a t depends on m only. This is Dr. Hill's primary solution or intermediate 
orbit. It contains only two arbitrary constants, while the general solution of the 
equation 9 (ii) contains four. It is desired to know the general solution of this 
equation and also of equations 9 , 9 (i) , 9 (iii) . We shall for simplicity at present 
take the equation 9(i), the extension so as to include the terms dependent on z 
being perfectly evident. We may write it 

xu =-Bu-Gs, 



(D + mfu + im 2 u + f»i ! s- t^ = 



where the terms on the right-hand side are of the order e' at the lowest. If now 

we put 

u = u + 8u , 

s —s -\-8s , 
u , s satisfy the equation 

(D + mfu -f im\ + t»» 8 *o-(Sl = °> 

and 8u, 8s will therefore be determined by means of the equation 

(D + my 8u + $m*8u + |m'fc - [ ( - + fa) f^ + fe)i - ^] 

= — B(u + 8u) — C{s + 8s). (12) 

We now suppose that 8u, 8s are small enough for the expression in 
square brackets to be expanded in powers and products of these quantities. 
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The expansion of this portion is 

1^ hu 3 JSs t _3_ (hif ^ _16_ (5s) 2 

8 s* 



x r i ou 6 os _£_ [puy 
u\sl L 2 m ~~ 2 s„ ' 8 «*o 

+ 4 « n sn 16 ul 16 s2 16 «gs 16 u Q sl + J 



As u , s are supposed known, we can, from their values, obtain the coefficients of 
the various powers of 5m, 8s in the form of known series. It is understood that 
the numerical value of in is always used, so that these series can all be expressed 
as odd or even power series in £ with numerical coefficients. Since 

u ^~ 1 = a Xa i ^ i , 
s £ =a 2a_ i £ 2i , 

all the series can be expressed in terms of even positive and negative powers of 
£ , so that, by a suitable arrangement, the largest coefficient is that of £°. It is 
assumed, unless otherwise stated, that the summations include all positive and 
negative integral values oft, including zero. 
Let 

x b — yp /a x s — vp/2» x s — ^nW 1 

«s»o nsi «m 

x b — yp«> x s — ypya 
—5—5 — *-"<s ' — r~T — "*-*<s i 



~^=:w , 4-i 2 =^-i^. 






When £ -1 is put for £, u changes to s and s to w . Hence Jf i =iC i and 

As for the calculation of these, M t is immediately obtained from the value 
of x/rl given by Dr. Hill (I, p. 249). For N t we have 

xZT % —*_ ** a £~ 8 

1 5 "■' o • o • 

Mssg rg r% 

"We calculate 1/r 2 and the real and imaginary parts of w 8 £~ 2 . These are easily 
done by the method of special values, which method will also give quickly the 
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required coefficients. In the same way we may find the other coefficients 
P, Q, R . . . . The whole process does not entail very much labor. 
The following table gives the values of M , N, P, Q: 



M 


= + .58902 


22856 


4 


M ±3 = -t- .00000 


06029 


7 


M ±i 


= + .00630 


84231 


2 


M ±i = + .00000 


00056 


5 


M ± > 


= + .00006 


28883 


4 


M ± 5 = + .00000 


00000 


5 






tf = + 1.75707 


88032 7 






Ni 


= + .03686 


55171 


8 


#_i = + .01078 


63527 


2 


^ 


= + .00054 


79401 


6 


N_ 2 = +.00001 


25690 


4 


N, 


= + .00000 


70129 


7 


N_ s = + .00000 


00982 


3 


N* 


= + .00000 


00824 


6 


JWL 4 = + .00000 


00007 


6 


N* 


= + .00000 


00009 
Po = 



: 1.17156 


77322 






Pi 


= + .02280 


40093 




P_ 1= + .01084 


18484 




P» 


= + .00032 


38766 




P_ 2 = + .00010 


24640 




P 3 


= + .00000 


40164 




P_ 3 = + .00000 


09526 




P* 


= + .00000 


00465 




P_ 4 = + .00000 


00092 




P* 


= + .00000 


00005 

$0 = 


= 1.17132 


P_ 5 = + .00000 
34260 


00001 




Qi 


= + .03476 


15314 




Q_x=— .00112 


12092 




Qz 


= + .00066 


73632 




Q-. z = + .00000 


31923 




Q 3 


= +.00001 


04714 




Q_ 3 = + .00000 


00337 




Qi 


= + .00000 


01462 




Q_i = + .00000 


00007 




Q* 


= •+ .00000 


00020 











The numerical values of these coefficients have been given here because 
they are fundamental for this method of development. Further, they are the 
same for every series of inequalities which it may be desired to calculate. The 
values of M, N, P, Q given above will suffice to determine all terms in u, s, 
depending on e, e 1 , I /a', y 2 as far as the first and second powers and the products, 
two at a time, of these quantities are concerned, and in z of the same classes of 
inequalities each multiplied by y. These form by far the larger portion of the 
expressions given by Delaunay. 

When we wish to determine the terms dependent on a' there are certain 
parts of the disturbing function to be added to the right-hand sides of the equa- 
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tions. But as these are homogeneous functions of u , s of the second, third, .... 
orders, their expansion is quite easy, and in order to keep the exposition as 
simple as possible they have not been written down. 

The various cases will now be considered in the order which seems most 
suitable to the objects in view, namely, the development of the theory and the 
special theorems to be treated in Part II. 

3. The Terms whose Coefficients depend on m, e' only. 

We take equation (12) and it is desired to find a particular solution of it 
such that the coefficients depend only on m , e 1 . As the primary solution depends 
only on m, and as the introduction of ef brings with it multiples of the angle 
n't + s\ no new constant of integration will be introduced. Also ou, Ss will here 
contain e' as a factor. 

At the beginning we neglect powers of e', and therefore of 8u, Ss above the 

first, and put for these terms 

o~u = u n , #« = $,,, 

and the equation may be written 

The coefficients B, G contain t only in the form ^(n't + «') or in the form £ ±m . 
(The constant part of the angle may be omitted since we suppose that m is not 
commensurable with any whole number, and e 1 plays no part in finding the 
coefficients of the terms in «,,, s v .) The solution is therefore of the form 

where j takes all positive and negative integral values including zero. 
The equations which determine the coefficients m > n'-j are 

(2/ +1 + 2mf n^Mm-i + StNtyA-j = coeffi of ?i+* in -2, [Ba } + Ga^^, 
(2j - If nLj +Stf-dui +X i N_ iVj _ i =: coeff. of r 2j '- m in -2, [5a,+ Co^J ?. 

From these we can find ty, qLj in terms of m by continued approximation. 
The next step is to find the terms of order e 1 ' by putting 

8u = u r) + u 7? , Ss = s n -\-s n i> 

(where the terms in u, s of the order e'' are u n ,, s^), and separating out those 
42 
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terms which are of the order e'\ The equation for them becomes, according to 
the previous notation, 

- coeff. of e' 2 in - ^ [B (u + «,) + (s + «,)] . 

The terms on the right-hand side and the known ones on the left-hand side 
necessarily contain £ in the form 

the solution is therefore of the form 

tv = «o«j C(w')i P + (w)i £ 2i+2m + W)j ^'- 2ro ] , 

and a similar expression for s^ . The coefficients (5727% etc., are the numerical 
quantities to be found. The method of procedure is the same as before. There 
are certain multiplications to be performed in the calculation of such a term as 
i£""VS SP-i^ 84 , but they can be performed in two steps by first computing 
^u 7t .'ZP-i^ ii (this result being wanted in the determination of the terms of order 
ee') and then multiplying the result by £w, . 

We can thus proceed to find the terms dependent on all powers of e', m by 
successive approximation and arrive finally at a solution 

u e > = ^o + m, + %> + • • • • , 

U e i := S -f- 8^ -\- S v t t • • • • 

It may be noted that the notation used here is intended to be a suggestive 
one and to point at once to the order of the terms and to the power of £ to 
which they belong. The letter 37 is associated with e' or the index m of £. In 
the coefficients, q, is the coefficient of £»•+» in u, (9797)^ or (^ that of £*+»", etc. 
This method of notation for the coefficients is found to be a convenient one and 
it is used for the inequalities all through. For example, if 2g be an even positive 
number, 

+ {n'% ¥ s+2qw + (*A £*-*"] . 

In particular we must notice that the terms in u which do not contain the 
angle n't + e' in their argument are 

a,S, [a, + e" (W)i + e' 4 W)i + ...-] £"• (13) 
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We can suppose now, for theoretical purposes, that all the terms dependent on 
e 1 , m have been calculated with the necessary accuracy, and by using the numeri- 
cal value of e' in (13) we can take that forjfche primary value of w instead of the 
value (11). But in the actual developments of the lunar theory this simplifica- 
tion has certain disadvantages. It is easier to find the extra small terms due to, 
for instance, the difference between (13) and (11), than to calculate the corrections 
to M , N which would be necessary. Hence we commence, when finding the terms 
depending on the lunar eccentricity or latitude or the parallax of the Sun, by 
using the primary orbit defined by the equations (11). 

5. The Terms whose Coefficients depend on m and the Eccentricity of the Moon's 

Orbit. Motion of the Perigee. 

In accordance with the remarks just made, we omit the terms depending on 
e', and in those results which are of a general character we shall merely state 
how they may be included and in what way the various expressions are altered 
by their presence. 

We take the equations 9 (ii) and putting 

u = u + u e + u e , + . . . . , 

s = *o T s e ~T s e* "T • • • • » 

where u , s have the values (11), obtain for the determination of u e , s e 

S~ l (D + mfu e + f- 1 u^M^ + fo.'SNt? 1 = 0. 
The solution of this equation is known to be of the form 

f- 1 ^ = 0,3 [%£»+• +«#»- ], 

Since c is supposed incommensurable with any whole number, the constant of 
the argument, an arbitrary of the solution, is not put into evidence. It can be 
restored at any time if necessary. This remark applies generally. 

When the substitution is made in the equation of motion, we obtain a set of 
linear equations to find Sj, ej without constant terms, and this fact necessitates 
a relation between them which determines the quantity c. Since the extension 
of Adams' theorems proved in Part II, and the results in Part III, depend on 
this relation to a large extent, it is necessary to develope it somewhat fully and 
to see in what way c is involved. 
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(14) 



Equating to zero the coefficients of % 2 ' ±e , we obtain the series of linear 
homogeneous equations in s j} e-, 

(2/ + 1 + m + of % + XiMtSj-i + 2M-j = 0,1 
(2/ _ 1 _ m + c )» e i. + ZiM^eLj + 2, #_,*,_,= . J 

The determinant formed by the elimination of e,-, e- we may denote by 



A(o) 



where the cross lines denote that d, d(c), d( — c) merely occupy the positions 
in the determinant assigned above. We have put d(c) for 



d(c) 


d 


d 


d(—c) 



I (5+m+c) 2 +J/ , M_ x 


M-> 


jc, 




if., " , : 


M x , (3+m-fc) s +if 


M_ x 


if_, 




M_ t , : 


if, , if 


(H-m+c) 2 +Jf , 


if-! 




ml* , : 


: Jf. jr. , 


M, , (- 


-l+m+ C ) S +if , 




id 


• 


a; 


*i 


(- 


-3+W2+C) 2 +if» , * 



and eZ ( — c) for the same expression with the sign of c changed. Also d stands 
for 



: Jsr , 


tf-i, 


M„ 


#-„ 


tf-4, : 


: -^i, 


^ , 


iLi, 


#-•, 


Ms, • 


: #„ 


^i , 


iVo , 


^-i, 


#-„ : 


1 ^8, 


^ , 


■tfi , 


#0 , 


«.i, : 


: N it 


■N. , 


^ , 


Xl , 


iv„ , : 



When these are placed in the expression for A (c) in the positions assigned and 
the cross lines taken away, we have the required determinant. 

This is a third form of the infinite determinant giving c. It differs from 
that obtained by Dr. Hill (II) in being what may be called "doubly-infinite," 
that is, it is infinite towards the directions not only of the outer sides of the 
square but also towards those of the cross lines. In this feature it resembles 
that which I obtained (VIII) by the use of the equations (8), but differs in one 
most important aspect, namely, that by a suitable system of divisors it can be 
put into a convergent form. 
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I repeat the observations there made as to the number of roots and the rela- 
tions which hold between them. Denoting by q the number of coefficients of the 
series e lt e 2 , . . . . e M , there are Aq + 2 unknowns in the 4^+2 equations; the 
determinant has Aq -f- 2 rows and columns, and it forms an equation for c of the 
order 2 (Aq + 2). If c lt c % be two roots of this equation not differing by an even 
integer, the whole series of roots is 

ioi, ±(c x ±2), ±(^±4), , 

± c 2 , ± (c 2 ± 2) , ± (c;j ± 4) 

Either c x or c 2 is shown to be zero, since 

u e = Du , s e = Ds 

is a solution of the equations for u e , s e . 

Now these two sets of roots are also roots of the equation 

( S in* J". - sin 8 ^)(sin* ^ - sin 2 ^) = , 

or of the equation 

(cos nc — cos TtCjXcos nc — 1) = . 
Hence we have 

A (c) == A (cos nc — cos nc^ (cos no — 1), 

an identical equation in which A is independent of c. We most easily find the 
constant A by equating the coefficients of the highest power of c in each member 
of the equation. In A (c) this is unity, and in the right-hand side it is 



4 * 4 s 4 * 

1? ' ~# ' 2Y 2 



/4 4 4 \ 4 

VT'375 * O ""J 

and therefore 

(cos nc — cos ?tc 1 )(cos nc — 1) = A (c) X f— . 8 • 8 ... .J 

=V(c). 



The new determinant v ( c ) i g P u t m to a convenient form by multiplying the 
middle rows of the upper half of A (c) by 4 , the two rows on either side of this 

by —, , the next two rows above and below by and so on. The lower 

4 — 1 8 — 1 

half being treated in the same way, all the divisors are accounted for. "We might 

then proceed as Dr. Hill has done to obtain the determinant in a better form, 

with every constituent in the central diagonal unity. The new determinant and 
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V (c) are both convergent, that is, in each case the product of the constituents of 
the central diagonal is finite and the sum of the other constituents is finite. (Poin- 
car6 IX, Helge von Koch X). But since the determinant will not be used for 
finding c, and since it has been treated as far as is necessary for the purposes of 
this paper, it does not seem to be of advantage to develope it further. Two 
points are of importance in what follows: firstly, the determinant A(c) is a 
symmetrical one, and secondly, in consequence of this and also since it van- 
ishes for the special value of c which denotes the principal part of the motion 
of the lunar Perigee, the coefficients e,-, ej are proportional to the first minors of the 
constituents of any row or column. 

It is not difficult to see what changes are made by the introduction of the 
terms dependent on e" and its powers. Instead of using the values (1 1) of w , s , 
we must use the values (13) ; also, instead of M t , N t we shall have 



= M ( -{- terms containing d* and its powers, 

(N) i = N i + " 

No other alteration is necessary in the equations since the terms arising from 
A, B, G [equation (3)] can be included in (M), (N). We, however, do not 
have (M)t = (M)_i. All the results concerning the infinite determinant for c 
still hold. Let for brevity 

(e) 3 = ej (1 + terms containing e" and its powers) , 

(«fy = *i(l+ " " " " " ). 

(c) = c + terms containing e'" and its powers 

be the new values. Then as before (e),-, (e'),- are proportional to the minors of 
the constituents of any row or column of the new determinant, which remains 
symmetrical. 

It remains to state that when c has been found, all the coefficients e t , e[ can 
be found in terms of one of them by any suitable method of continued approxi- 
mation. The minors of the determinant do not seem adapted to calculation. It 
is found that e — s' is approximately 2e where e is Delaunay's eccentricity, and 
therefore this quantity is a suitable arbitrary. The letter e is used in all coeffi- 
cients involving the lunar eccentricity just as the letter ri is in all coefficients 
involving the solar eccentricity. 
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5. Further Approximations to the Terms dependent on the Eccentricity of the 
Moon's orbit and to the Motion of the Perigee. 

The next process is to find the terms of the order e 8 . These follow in the 
same way from the results just obtained as those of order d* did from those of 
order d. If the part of c depending on e, before neglected, be he, we can easily 
see that he is of the order e a and that we need not know it to obtain the coeffi- 
cients of the order e 2 . We suppose the coefficients of the order e 3 found. 

We proceed, then, to find the terms dependent on e 3 with the corresponding 
new part of c. We are only going to consider the terms depending on the 
indices 2/± (c + he) of £, since it is by these terms alone that he is determined. 
Stopping at the third order we put 

u = u + u e + u e , + u e3 , 
where 

u e + u e , =a x J [^ +e+ic + <?;£*—* + ( 6 v) i ^+ c+Sc + («v- to ]. 

and similar expressions for the complementary variable. The last two coeffi- 
cients are of the order e 8 , and therefore on substitution in the equation of motion 
we can put c for c -f- he when this quantity appears as a factor of these coeffi- 
cients. We obtain, by separating out the terms of the order <s s , 

[(2/ + 1 + m + o + hef - (2/ + 1 + m + cf\ Sj 

+ (2/ + 1 + rn + of (<=¥), + X t M ( (*¥),_< + 2 t N { («=")*_, 

= coeff. of ^+°+ s ° in ™L— (u + u 6 + ^) _i ( s o + s e + s 6 ,) -i of the order e 3 , 

a$ 

and a similar equation. Omitting negligible quantities on the left-hand side, 
the equation becomes 

28c(2j^l + m + c)e j ^(2j + l + m-\-cf(^) j +X i M i (^)^ i +^ i N i (eni-j 

= same expression as before. (15) 
Similarly the other equation gives 

2he(2j— i- m +c) E >_j+(2j-l— m+cXen-j+ZiM-deOi-j+ZiN-ii^j-i 

= coeff. of ^-v-— * in ^ (w + u e + w e2 ) _i (*o+ #e+ Se») -1 of the order e 8 . (15) 



a, 



o 



The terms on the right-hand sides of these equations are known, those on the 
left-hand sides are to be determined. We can get an equation of the first order 
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for Sc free from the as yet undetermined coefficients («V), (ee' 2 ). There are, in 
the notation of the last section, 4q + 2 equations and 4q -f- 3 unknowns to be 
determined. One of these is arbitrary, and it involves a new definition of the 
"eccentricity constant" of integration. We may define this as may be most 
convenient. At present we proceed to the determination of Sc , which is not 
affected by this observation. 

Multiply the first of equations (15) by e jt the second by sL j} and sum the 
two for all values of/. Now, in these equations the coefficients of (e 2 e% (ee 1 ')^ 
are exactly the same, term for term, as those of ej, eLj in equations (14). If, 
therefore, we multiply the equations (15) by the respective minors of A(c), all 
terms involving (eV) j( (ee' 5 )-j will disappear. But as these minors are proportional 
to e f , sLj, the same object will be attained if we multiply by e f , eLj. If, however, 
we conceive of difficulties arising from the fact that A (c) is not in a convergent 
form and therefore that its minors are infinite quantities, though their ratios are 
finite, we can proceed as follows. After the multiplication by sj, sL, and the 
addition of the equations for all values of j, the coefficient of (eV)* is 

(2» + 1 + m + cf Si + XjMje^ + 2,^ef_ <f 

which, by equations (14), is aero. And so for all the other similar terms. This 
process is, of course, algebraically identical with the previous one. 
Hence we have 

2&2, [(2/ + 1 + m + c) 6 ? + (2/ — 1 - m + c) <,] 

= 2, fs, X coeff. of £*+« + "• in &2 (u + u e + VT* (s + s e + «„)-' 

+ bLj X coeff. of £-»—» in ^— (u + u e + «„)-* (s + * e + «•.)! , 

where the terms on the right-hand side of the order e 4 alone are considered. -But 

and therefore, using this on the right-hand side of the previous equation, we 
obtain 

2&?2, [(2/ + 1 + m + e) e? + (2/ - 1 — m + c) £,] 

= const, part, order ,, in £ . S + 4VC& + ?+.,)■ ' (16) 
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Transforming to real coordinates and putting 

X* — xo + x e + x e i , Y 6 * — y Q + y e + y* , R% = X% + 75. 
the right-hand side becomes 

const, part, order e*, in expansion of — g • ™ — , (.lo'J 

an expression of somewhat remarkable form and simplicity. When first ob- 
tained, it seemed to suggest that Adams' theorems (IV) as to the connection 
between the constant parts of the Parallax of the Moon and certain parts of the 
motions of the Perigee and Node, must really arise naturally from this mode of 
development of the lunar theory. After some trouble I succeeded in proving 
again these theorems and extending them so as to show them not merely in the 
form of ratios but in an exact form. They will be given in Part II of this paper. 
It is easy to extend these results so as to include all powers of e 1 . "When 
we do so the formula (16) or (16') becomes 

23(c)2,[|2y+ 1 + m + (c)} (e)f + (2; - 1 - m + c)( e ')LJ 

= const, part, order e\ in \ . ( x )e^l + (Y) e ,(x) e ^ (1?) 

where the brackets round the various symbols have the same meaning as that 
given at the end of section 4. 

The coefficient of ho or h (c) remains the same if we are finding the terms in 
the motion of Perigee dependent on the latitude of the Moon or the parallax of 
the Sun. If we are in the process of finding h (c) , that part that depends on e' 3 
can be obtained as in equation (16) by finding the constant of order eV 2 in the 

expansion of 

J^_ (X\x, + (Y) e y e 
at ■ (R)% 

where (X) e , (T) e , (B) e contain merely the terms, supposed to be previously 
found, dependent one, ef, d\ ee' and m. There are also the terms arising from 
6) 2 to be included. So that in all cases the determination of the new part of c 
depends only on an equation which contains known terms found in previous 
approximations. 

Similar results will follow for the terms dependent on higher powers of 
e 2 , e' 2 in c. For example, the equation necessary to find the part in c of the 
43 
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order e 4 may be easily obtained by looking at equations (15). It is <5 2 c where <$ 2 c 
is given by 

WcXj [(2/ + 1 + m + c) $ + (2/- 1 — m + c) <£,] + (Sc) 8 ^ [ef + el',] ^ 
+ 2 j [(2i+l + m + c) 2 ej ( £ V), +(2i-l-m + c) 8 ^,( e6 ' s )_J I (lg) 

= const, coeft., order e 4 , in — ^ • ™ > 

a xt e 4 J 

where X ei = aj + x e + ic e2 + a5 e> + a e 4 , etc. 

The coefficients to be determined therewith by the equations from which this 
expression arises, are those of the order e 5 with arguments (2iD± I), where 2D 
is the argument of the Variation and I that of the principal elliptic term. 

If we wish to determine the terms dependent on the square of the parallax 
of the Sun, similar results will follow. On the right-hand sides of the equations 
we shall have further terms arising from IX X , but these terms will be all known 
by previous computations, so that the general result for all terms in the motion 
of the perigee is this : In any stage of the approximations we can always find the 
new part of the motion of the perigee by simple computations without being compelled 
either to solve an infinite determinant or to find by successive approximation the new 
coefficients arising therewith. The new coefficients can then be found by the usual 
methods. 

These results hold for any terms whatever in the motion of the perigee, and 
also, as will be seen immediately, for those in the motion of the node. The 
most important feature in them all is that just stated in italics. The solu- 
tion of an infinite determinant, in whatever way performed, is very laborious, 
and the continued approximation method applied to &c and to the coefficients 
simultaneously, is not much less so. The method here given, when once the parts 
dependent on m only have been calculated, avoid these laborious processes — an 
important point quite independent of its mathematical interest. It is not to be 
forgotten also that in one respect theory here goes hand in hand with observa- 
tion. The motions of perigee and node are now capable of being determined by 
observation with much greater accuracy than the majoiity of the coefficients of 
the periodic terms. In the method outlined above theory also gives the new 
parts of the motions before the coefficients of the corresponding periodic terms. 

We now proceed to find the effect of including in our equations the motion 
of the Moon in latitude, 
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6. Terms dependent on the Latitude of the Moon. Motion of the Node. 

For these terms we use equations 9 (iii), omitting those dependent on the 
solar eccentricity. Since z is a small quantity in powers of which we 
suppose that expressions are possible, and since the second of these equations 
contains it in its first, third, fifth, .... powers only, while the first equation 
contains it in its second, fourth, .... powers only, we are able to treat the 
equations separately, going first to one and then to the other as we obtain 
each new approximation. This fact facilitates the performance of the actual 
calculations very largely. 

Let y denote generally a quantity of the same order as z. The s-equa- 
tion is 

Dh — (m* + * J ) 8 = 0. 

\ (f + zyJ 

Let z = z y + *V + V + ••••>") 

u = u + u y »-\-Uyi'+ > (19) 

s = «o + *y> + V + ' 

Expanding the last term in the ^-equation, we have 

n » / , i * N \ 3 xz 8 . 15 rcz 5 . 

Irz — [m? -\ — 5-3 = — -r + — r + •••• 

\ r 3 / 2 r 5 8 r 

The equation for z y will be 

D%-(V+^-) Zv = 0. (20) 

The solution of this equation can be put into the form 

ay=T=aoS [K^+o + KLj^-o-], y= + »....- «, (21) 

where g must be found, and K jt KLj are of the order y and such that K 5 = — KL S . 
The arbitrary constant attached to the new angle which the presence of g intro- 
duces can be understood, as in section 4, to be present. 

It is known from Adams' work (III) that g is found from an infinite deter- 
minant. Expanding x\r\, we have from section 2 

D\-2z y 2M i ? i =0, 

and therefore on substitution of the assumed value of z y and on equating the 
coefficients of ^ 2i+ " to zero, 

(2j + gfK j — 2XK i M j _ i = 0. (22) 
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This set of equations gives a symmetrical infinite determinant which, after 
certain transformations, was calculated by Adams (III) to find g. The coefficients 
Kj are proportional to the first minors of the constituents of any row or column 
of this determinant. When g is known and we have thence obtained the ratios 
of Kj to any one of them, say K , which we take as the arbitrary constant of the 
solution, we substitute the resulting values of z in the first of equations 9 (iii) and 
find by the same methods as before the values of u yi , s y *, which, as the suffixes 
indicate, are of the order y % . 

Returning with these to the z-equation, it is desired to find the terms in z 
of the order y 3 together with those in g of the order y % . Confining our attention 
to the latter object, we assume 



"We omit the terms which depend on the indices 2i± Zg, as their coefficients 

can be determined to the order y 3 with the known value of g. We have 

(JPK')j = — {KK")^. Equating to zero the coefficients of ^+s+ s s for all values 

of j, we get a set of equations for the determination of the new coefficients and 

Sg . Let 

X y * = x + x y * , Y y , = y + y yi , E^ = X^ -f 7^ , (23) 

then, by reasoning similar to that used in the determination of Sc, the set of 
equations may be written 

2$g (2/ + g) K 3 + (2/ + gf {JPK%- 2XMj_ { {IPX'), 

= coeff., order y\ of ?+>+« in *%f~* . (24) 

Multiplying this equation by X jt and summing the set of equations so treated for 
all values of j, we have by equation (22) (as in the last section) the new coeffi- 
cients disappearing and 

28gX (2/ + g) Kf = [2 K, X coeff., order y\ of e + " + Sff in *5^=jL ] . (25) 
Similarly, 



2302 (2/ + g) KHj = S [KLj X coeff., order y\ of ^-,-H in _^L_A_] . (25) 
Now, since K'_ 5 =- — K jt and therefore 
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we obtain, by addition of the two previous equations, 

4Sg X(2j+g)iq= const, coeff, order y\ in „^L— , ( 2 6) 

which is the formula for obtaining Sg. 

This formula is easily extended so as to include all powers of e'. In deter- 
mining K jt instead of using r we use (r ), where the notation is that used in 
section 5. Let the corresponding values of K j} g, Sg, etc., be (K) jt (g), S(g), 
etc. The equation then becomes 

4o (g) 2{2j + (g) \ {Kf 5 = const, coeff., order y\ in fl , | ( ^ 8 ^. (g) ^, • ( 2 ?) 

The part of the motion of the node which is of the order e 2 , and that of the 
motion of the perigee which is of the order y % , will be required for the investi- 
gations in Part II. They are not difficult to write down if we look at the 
results (16') and (26). Let z ey be the part of z of the order ey. This can be 
found when we know the principal parts of the motions of perigee and node only. 
It can then be easily seen that the part of Sg depending on e 3 is given by means 
of the equation 

4hg Xj (2j + g)K?= const coeff., order e 8 y 2 , in ^vfo + O . ( 2 g) 

Also that the part of Sc of the order y % is given by the equation 
2&J2, [(2/ + 1 -f in + c) «» + (2/ — 1 — m + c) e'l^ 

— const nart order eV in - (X e + x yi )x e + (Y e -\- y y ,)y e , * 

_ const, pait., order ey,m^. [(J; - ^ )2 + (y . - ^ - ^^grji ■ W 

Finally, by putting brackets ( ) round each letter, we may include as before 
in the formulas (28) and (29), all terms dependent on d* and its powers. 

The equations (28) and (29), with powers of e" included, will be named (30) 
and (31) respectively. (30), (31) 

It is unnecessary to go further with these expansions. Sufficient has been 
said to indicate the manner of treating the whole- Lunar Theory after this 
method. The terms which depend on the parallax of the Sun are a little trouble- 
some owing to the new terms which arise from £i lt but they are much less so 
than when treated by the simultaneous equations (8) , since, for a given degree of 
approximation, they are one order lower in u, s in the former case than in the 
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latter. Nevertheless, should it be considered that in finding the terms of high 
orders in e, e 1 , y, 1/a! the method indicated above becomes too troublesome, 
nothing prevents us from returning to equations (8), these being available at any 
stage. No increase of labor results from the mere change of method. 

After the experience gained in the calculations of the coefficients made up 
to the present time by both methods, it seems probable that the equations (8) 
will involve less labor than (9) for terms of the fourth order in e, e', y % , 1/a', if 
not for the third. Nothing even then stands in the way of the determination of 
the new parts of c and g by the methods just exhibited. 

7. Order of Proceeding. 

Having resolved the chief difficulties in the theory, namely, the determina- 
tion of the motions of Perigee and Node, it is possible to give a consistent plan 
for the computation of the complete expressions of the coordinates. The order 
which seems to arise most naturally from the previous results may be sketched 
out as follows : 

Terms dependent on 

(i) m only. 

(ii) m and ef, e' 2 , e' 3 .... in succession. 

(iii) vn and e, eef, ee 1 ', . . . . , e 2 , eV, eV, . . . . , e 3 , eV, .... 

(iv) m and y, ye 1 , ye", . . . . , y*. yV , y 3 , . . . . 

(v) m. and any other combinations of y, e, e 1 . 
(vi) rn and 1/a' and combinations with y, e, e'. 

This order may be varied somewhat. 

The arguments in Delaunay's notation are, 2D that of the variation, I that 
of the principal elliptic inequality, I' that of the annual equation, i^that of the 
principal term in latitude. The expressions of the coordinates should then be 
arranged in the order 

2iD, 2iD ± I', 2iD ± 2?',- , 2iD ± I, 2iD ± I ± V, 2iD ± I ± 21', , 

and so on. The chief difference from Delaunay's arrangement being that the 
terms whose arguments differ only by multiples of 2D are to be placed together 
in all cases. 
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It is proposed to carry the computations sufficiently far in x, y, z to enable 
the coefficients in longitude and latitude to be found accurately to one thousandth 
of a second of arc. This makes it necessary to carry some of the earlier calcu- 
lations to a somewhat greater degree of accuracy than is actually needed for the 
special coefficients. But in any case the terms in in only, and the principal parts 
of the motions of perigee and node, carried to 15 places of decimals by Hill and 
Adams (I, II, III), are amply sufficient. 

8. The Arbitrary Constants. 

The constants which arise in the arguments by the integrations are the 
same as those of Delaunay. Those outside of the signs, sine and cosine have 
been referred to generally as a , e, y. As for a , we can define it in such a way 
that the determination of the coefficients is made as simple as possible, and this 
will cause in general a determination of a new part to a whenever we are 
making an approximation which contains a constant term. It would seem better, 
however, to define a as that value given by Dr. Hill (I) in terms of {fi/n % ) h and 
m , making the necessary changes after each approximation when the latter has 
been made under the previous method. Similarly, it seems the best to define 
the "eccentricity" constant so that it is the coefficient of sin? in x/a , and the 
latitude constant so that it is the coefficient of sin F in z/a . Then a differs 
from Delaunay's a by terms which vanish with m; the first approximation 
(neglecting m and all powers and products of the constants) to the eccentricity 
constant is 2e, where e is Delaunay's constant; and the similar first approxima- 
tion to the latitude constant is the quantity defined by Delaunay as y. 

The numerical value of m is used all through, the other constants being left 
arbitrary. 

9. Verification of the Results. 

For verification purposes we have several equations at our disposal. The 
most useful of these when the coefficients of %® are not under consideration, is 
the second of equations (8). If it be desired to verify the coefficients of ^'+ A , 
where A is some linear combination of m, c, g, we substitute in this equation 
the results obtained and find therefrom an equation of the form 

%£ A + */%-* = o, 
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where 2,, V contain integral powers of £ only. The coefficients of the various 
powers of £ ought to vanish identically after substitution, and therefore if we 
put £=1 in X, 7J, the results ought equally to be zero identically. When 
£ = 1 , /l + a/ = 0, and therefore this method will not be available when .A is 
zero. In this equation of verification all the terms are of the second degree in 
u, s, and the combinations of the coefficients are totally different from those 
arising from equations (4) used for finding them. 

As for the coefficients of £*■*, the results proved in Part II concerning the 
constant part of 1/r used in connection with equation (6) will verify all coeffi- 
cients of the orders e 2 , y 8 , e\ e 3 ^ 2 , y 4 and all powers e' 2 . For practical purposes, 
then, we have only the terms of orders e l% , e", . . . . lfa'\ e*/a'\ . . . . , etc., out- 
standing. Equation (6) can be used for these by expanding 2%{us + tf), and, 
after choosing out the terms of the required order, putting £ = 1 . 

10. General Observations. 

It will be well to contrast the method here outlined with that used in the 
development of the Elliptic and Parallactic Inequalities (VII, VIII). It must 
be premised that both methods are built up on Dr. Hill's primary solution. 
These remarks refer merely to the best method of continuation. 

The chief disadvantage of the present treatment lies in the frequent multi- 
plication of series of £ with numerical coefficients. Although every multiplica- 
tion is available for several classes of inequalities, besides the one for which it 
was especially made, in terms of high order these computations will become 
very tedious. But there are several advantages which may render the method 
worthy of consideration. In the first place there are no numerical coefficients to 
be computed after the multiplication of series has been performed, like those 
arising from the use of equations (8). In the latter, for instance, when finding 
the terms of order e 8 we have to compute 

2,0'. *'• °> l ) e i e l-i 

for all values of *', the portion in brackets being a somewhat complicated func- 
tion of m and c, and different for different values of i,j (VIII, p. 322). Such 
calculations are fruitful in producing errors. Secondly, the terms arising from 
£i x are very simple and, when I fa' is neglected, of the first order in u, s, so that 
only a very few short operations are necessary. Thirdly — perhaps the most impor- 
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tant point to be mentioned — the multiplications of series, being always of the 
same nature, can be performed by a computer from one general formula. This 
would require little extra knowledge on his part beyond that of logarithms, and 
his results can be checked by a computation with a special value of £. 

The method here outlined has been successfully employed in computing the 
coefficients of orders e 1 , ef, eel with the numerical value of m ;* much labor appar- 
ently was saved by the readiness with which the work was reduced to easily 
remembered operations and by the diminished opportunities of making errors in 
calculation. Certain difficulties owing to the near equality of the coefficients of 
£ ±m and ^±( a_m > and also to the large divisors, were more readily bridged. I hope 
shortly to test the terms in e 3 , and also the part of the motion of the perigee of 
the order e % (VIII, p. 337). This last differed somewhat from Delaunay's value 
(VI), and before it can be accepted, ought to have a further test applied and the 
errors, if any, corrected. 

PART II. 

Theorems concerning the Constant Part in the Expression of the 

Parallax op the Moon. 

In the investigations which follow, the Parallax of the Sun is neglected. 
With this exception the results obtained are perfectly general. 

Let x = X q , y=T g , z = Z q 

be the general solution of equations (9) or (10), neglecting quantities in the coeffi- 
cients beyond the order e p <y q ~ p (p = 0, 1, 2 . . . . q). We define e, y as the 
arbitrary constants of integration in eccentricity and latitude, of the same order 
of small quantities as those used by Delaunay. 

Let X q = oc -\- Xi -f- . . . . + x q , 

Y q= yo + yi + — +y q > 

Zq = Z l H~ ^2 +••••+ Zq , 
Rq = Xq + Y g + Af-1 > 

where x q , y q are of the form a q e q + a q _ i e q ~ !S y 2 + a (Z _ 4 e (2 ~ 4 y 4 +••••) 

and z q is of the form a g _ 1 e <z-1 y + a q _ s e q ~ 3 y 3 + . . . . 

From what we know of the expressions of the coordinates of the Moon, 
a g ,a a _2, .... are functions of m, e 1 and of the angles 2v(t — 1 ), n't + e', 

* Monthly Not. E. A. S., Vols. LIV, p. 471 ; LV, pp. 3-5. 
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(32) 
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cv (t — tj , 2gv (t — t z ) and a g _t, a q _ 3 . . . . of m , e 1 and the angles 2v (t — t ) , 
n't + e', cv (t — fj) and odd multiples of gv (t — t 2 ) . According to the notation 
previously used there should be brackets ( ) round all the symbols X, x, etc. 
These must be understood to be present, but are omitted here for the sake of 
simplicity of expression. We proceed to prove some theorems. 

11. If x, y, z have been calculated to the order 2q — 2, the constant part in 
the expansion of lj*/(x % + y z •+- z % ) can be obtained to the order 2q without further 
reference to the equations of motion. 

With the notation just outlined, this theorem may be more fully stated by 
means of the equation, 

3 



const, part, order 2q, in the expansion of 



R% 



= const, part, order 2q , in the expansion of 

[ 8 x X 2q _s + y Y H _, _ 3 g^, {xoXi + yoyi) j m 



(33) 



Jiiq — 2 



n 



The proof of this is based on Adams' method of proving his theorems referred to 
below (IV). The equations (10) are, if r 2 = a; 2 + y % , 



x — 2riy — Zn'*x + A'x + B'y = 

y + 2n'x 

z -\- n'*z + n' pz 



(IX 



(r*+z*y' 



(34) 



(IZ 



Let x , y , z be the values of x, y, z, expressed by the same formulae and 
with the same angles, but with different arbitrary constants of eccentricity and 
inclination (i. e. different values of e, y). We have then 



*o ~ 2n< y — Zn'\ + A'x, + B'y = - f?° »\t ' 

v*"o "T z o) 

y +2n'x +B ao +Oy 9 = -j^&^, 

+ ^ + .Vo =-(^)- 



(35) 
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Let the suffix correspond to zero values of e, y. We suppose as before 
that expansion may be made in powers of e, y. Since equations (34) hold 
generally, they must hold if, after the right-hand sides have been expanded, we 
stop at any given order in e or y in all terms. They will therefore hold if we 
stop at the order 2q in e, y. Hence we can put x = X 2q , y = Y iq , z — Z 2q _ 1 , as 
Zi q -i is the value of z correct to the order 2q. Multiply the first two of equa- 
tions (34) by x Q , y respectively and subtract the sum of the results from the 
sum of the first two of equations (35) multiplied by X 2g , Y iq respectively. We 
can put the resulting equation into the form 

-fa [X^o + Y iq y — x X 2q — y Y m — 2ri ( Y iq x — X n y )] 



= — ft (-X^aso + Y Zi y )(j^ ^ . 



Since the expressions of the coordinates contain t only under the signs, sine 
cosine, we deduce immediately, 
const, part to the order 2q inclusive, in the expansion of 

(X 2q x + F 2g2 / )(^ 7?) = ° • ( 36 ) 

Let the notation G iq [§] denote "the constant part, order 2q, in the expan- 
sion of Q " according to any parameter — which we can put into evidence when it 
is necessary for the sake of clearness. The equation (36) then gives 

q r X 2g a- + Y ig y -] _ G r X 2q x Q + Y iq y -i 
L M 2q -1 L Mq J 

_ q^ p^o + ^yo -] t ( 37) 

since the parts of X 2q , Y 2q of order 2q are x 2q , y iq and Rl=Xo+ Y$= xl+yl=. r\. 
The equation (37) holds for any positive integral value of q. 
Now, if we stop at the order 2q, 

R 2q — R^q-i + % x o x tq + 2y y 2q + 2z 1 z 2q _ 1 , 
and therefore 

— 3 {x\ + yfl "^ + ^ + ^3=l \ . (38) 
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Putting r% for x\ + y 2 , we deduce from (37) and (38) 

q r gojg izz jj+yojVri ~i _ q r 3 goggg + yoya? + gi%g-i "i # ( 39 \ 

Again, by expansion, 

r r i "I _ r> r L_ «ng» t + yoy«g + g|g^-i i 

2ff Li? — J ~~ 2q L7? — — r s J • 

Hence by (39), 

q r 3 "1 = (7 r 3 _ ^-^u-i + yo^tg-i i 



_ r r ? 

2g— a ~r yo-*22 

(n- 2 +2^o, a2 _ 1 +2^ )2ff _ 1 )5 



^0-^29-3 + yo^-2 + #0, 22-1 s l g gg— 2 ~j 



■ (33') 



where we have put 

So, 3g-l == X X 2q — 1 T" VoVm — 1 "T" Z l Z 2g-2> 

Si, 22-1 := ^l^a-i + yiy2g-i + %2 2(2 _ a . 

Hence by expansion we obtain 

rt r 3 "1 /7 T 3 I Q O , 1 A , 2 ? - 1 o ^1,22—1 ^0-^22 — 2 + y0^_22-2 

^Lii^J L -«'22_2 r r - n '22-2 

I oO "0, 2ff — 1 Z 1 S 22 — 2 

+ 3/V 1 -j 

'0 
/-S I ,0 ^ q "O, 22-1 I ic "0, 22-1 Sq, l q »1, 2<?-l \"| 

w+^0,1^ 3-^- + 15— — ^ 3 -^;j 

r> r 3 X JC 2q — 2 "H VoXiq—i o Z l z 2q—2 M), l M 

— °8s Ip ds d 3— J J ' 

which formula, since s = , proves the theorem. 

12. Adams 1 First Theorem. 

We can deduce this theorem immediately from the above result. It states 

that in the constant part of the Moon's Parallax all the terms of orders e l and y* are 

zero. 

In equation (33') let q = 1 . "We obtain, since z = , 
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First, let y = ; this shows that the term in e s is zero. Next, let e = , then 
similarly the term in y s = 0. Also it is to be remembered that e, y occur only 
in the const, part of 1/B in even powers. In the above, e 1 * and its powers are 
included. 

13. Adams' Second Theorem and its Extension. 
This theorem may be stated as follows : 

V 0i [JL-] = Ee* + 2Fey + Gy* (40) 

and O t [(c)] = He* + Ky\ 

C i i{gy\ = Me*+Ny\ 

where E, F, G, H, K, M, N are functions of m, e 1 only, then 

K — F ' N ~ G" 

The proof of this theorem may be found in Adams' paper already referred to 
(IV). It is based on the consideration of the orders of certain functions. In 
the extension of this theorem which follows, I shall prove again Adams' 
theorem by another method which, though lacking altogether the neatness of 
his demonstration, has the advantage of putting the connections between the 
constant part of the Parallax of the Moon and the motions of the Perigee and 
Node in a rather clearer light. 

Extension of Adams' Second Theorem. 
Let 

<?T e = 2&2,tW + 1 + « + (P)} ( s )*+\2j-l-ni + (c) }(*?-,] , 

fT y =±^ } {2j+{g)-\{K)) 

(see equations (17) and (27)), where T e , T Y are expressed in terms of the same 
eccentricity and latitude constants (which we have called e, y) as 1/B, (c), (g) ; then 

ET e =6E, MT y =6F,\ ( . 

KT e =QF, NT y =QG,\ y } 

and thence T t M 

T 1 -~K' 
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As before we omit the brackets which denote the presence of powers of e', for 
the sake of simplicity of expression : it is understood that they exist in what 
follows. 

We have from equations (17) and (31) in this notation, 

e*ET t = G et p^+y.SV ] f 

j^ kt — a r x » ( x ° + x ^ + v* i r < + gg) i 

^ e ~ eV Lu^ e + a^)" + (F. + </ v ,) 8 + («, + -grJ * 
These can be included in the equation 

+BT. + ttfX*. = «-. » [pftS • («) 

For since x, y contain only even powers of y , Xi = £c e , yi = y e > and when y = , 
X 2 == JT 6 ,, etc., it therefore includes the first equation. If we want the second 
equation we neglect powers of e in X, T above the first, hence X 2 = X 6 + <*V> 
etc., for this case. 

Again we have from equations (30) and (27) 

*y*M T y = 0++ [ z r( g v + g ^)] , 

and these can be likewise included in the equation 

eYMT y + yW v = <7 ev , 7 , [^] , (43) 

since z contains only odd powers of y . 

Now, since constant terms contain only even powers of e, y, we have by 
addition of (42) and (43) , 

+BT. + <?f (KT e + MT y ) + y*NT y = A [ ^+^+ft^ ] . (44) 
By means of this equation we can reduce the four equations (41) to two, namely, 

For, suppose (A) and (B) to be proved. In (A) put y = 0, then 2=0, and by 
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(44) and (40) the first of the relations (41) holds. Next, suppose e = 0, then 
in a similar manner the last of relations (41) holds. Finally, consider the terms 
in (A) of the order e 2 y 2 - Owing to (B), the second and third of equations (41) 
will hold. Hence all that is now necessary is to prove the two equations (A) 
and (B). Since the methods of proof may be of value in investigating similar 
relations for higher powers of e, y in the parallax than the fourth, they are given 
separate statements. 

14. To prove that 

n f 6 1 a V x^ + y^ + ZjZj ] (X , 

° 4 \-7tJ ~ c * L — (izj+a**)! J ' (A) 

Attach to x x , x 2 , y x , etc., a factor a with an index of the same order in e, y as 
the quantities to which they are attached. Thus for x x we put ax x , for y%, afy 2 , 
etc. ; a will be omitted where its presence is unnecessary. Let 

o rs — X r X s "T" y r ys' 
Then B\ = r% + a?S lx + a 4 £ 28 + 2a£ M + 2a 2 ^ + 2a 3 £ 13 + a 2 z\ . 

Hence, reducing to a common denominator, the expression 

_6 2 ^ + mY * = -Jr (2^ +3a 2 ^ 11 + 3a 2 ^ + 3a<%+ 5aS 01 + 5^^+ 6a s S n ) 

= -i- (4r z +5a?S 11 + 5ah\+ 6a 4 £ 22 + 9a£ 01 +10a 2 £ 02 -flla 3 ,Si 2 ) 

5 4a 7 ^+5a 1, /y il +5aV 1 +6a 11 ^+9a 8 /y ol +10a !, /S f oa +lla 10 AS , i a 
(aV + «%+ aN + a 13 # 22 + 2a 9 £ 01 + 2a w S 0S + 2a n S Xi f 

xiXt + yiYs + 4 



a 



+ 



m 



— „» d ( 1 \ x x X % + y x Y % + z\ 
~ a da KBaO^ B\ ' 

Now, owing to the definition of a we can consider 1/B 2 as expansible in 
powers of a, and it is easily seen that 

d ( 1 



da \B 9 a 4 / 
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contains no term of the fourth order in a. Hence 

Gi 1 4- - 2 g^+jMSl = a 4 fgiA±|^±i] . (45) 

-O^ Xtjj - 1 Xl/2 

Let now # = 2 in the theorem of Art. 11. Equation (33) becomes 

° 4 L^J — ° 4 L^- rs 3 -^--l- 

Hence from (45), 

° 4 LsrJ - °* L ^ 6 -jr-J 

-°*L (22- +*»,%)• J» 
which is the result required. 

The proof of equation (B) will be divided into two parts, of which the first 
contains a method of general application to the terms which involve both the 
constants e, y. 

15. Lemma. To prove that 

G rx y *X e > + y v ,F e , _ x e ,X y , + y e ,Y y ,l _ Q ^ ,^ 

L ±t e -X My* J 

Let, for the purposes of this proof only, 

£ = % + a> + ««», £' = «b + av , 

t = yo + y e + */«»> *t = 2/0 + sv 

r 8 = f + »7 2 , £'=* v 

Then £, ^ contain the terms in x, y, independent of y, as far as the order e 2 , and 
£'> Yl', £'> ^i the terms in x, y, z, r, independent of e, as far as the order y a . 

The proof of the lemma is not difficult. We proceed as in section 11. 
Taking the equations of motion (34) and substituting therein for x, y, z first 
£ , ri , and then £', yf, £', we have 



n+ln'Z +_B>Z+C> V j = -n« r , 
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and e— W- 3nT+ A'j? + B' n '= — ii-fjr, 

#+ 2n>? + B$ + Ort= -p-£r', 

the right-hand sides of the first pair being expanded as far as the order e 2 and 
those of the second pair as far as the order <y z . Take the sum of the first pair 
multiplied respectively by £', vf from the sum of the second pair multiplied 
respectively by £, vi. The terms on the left-hand side of the resulting equation 
form a perfect differential, and therefore contain no constant part. Hence those 
on the right-hand side as far as the orders e z , y 2 contain no constant part. We 
have then in particular 

Since o contains no terms in j/ 2 and o' contains no terms in e 2 , we have 

°«v L — p — 7 s — J — u > 

whence follows immediately the equation (46). 

From this result we can deduce the result (B). This seems to be most 
quickly proved by simple expansion. 

16. To prove that 

p P^L+llZ?! — n P> Z *\ 

Let S el = x x e + y y e , S^ = x x y , + y y y , , 

S m = x x e ,+ y y e , , S 12 , = x e x y , + y 6 y y * , 
$n = x l +yl i S iV = x ei x yl -)ry ei y yl . 

The orders of any terms will then be denoted by the suffixes of S. We have, 
with this notation, to the order required, 

B y , = rl+2S 02 , + 4, 
and therefore, by expansion of equation (46), 
rt r$sz' o $oi $12' i ,p f o $02 + h $ii i 15 $oi 1 



rl 



45 
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or, after reduction and multiplication by 2, 

O^, [— 6 So1 f »' — 3 So2 '^ u + 15 ^ f oy + 3 ^^] = 0. 

*o *"o *0 ^0 

Rearranging, this equation may be written 
0*. [S 01 j - 3 -^ + 15 S oi^+iz\S 01 ) 

- 3 ^ii — * *^ ~ 3^y A- ^ i _ 3 ^ ±Mi + JJ- J^-1]=0, 

which is the expansion of 

^ T tin. + ^ u + JS& — A "I _ 

°«V Ltf + S n + 2# 01 + 2£ 02 + 2S W + 2S 12 , + z\p • 

Referring to the notation defined by equations (32), this is evidently the 
same, if we omit superfluous terms, as 

r r ^i^a + yi^a — gf ] — n 

°eV L ^3 J — u ' 

This equation is easily shown to be the same as 
whence follows the equation (B). 



The above results give some insight into the nature of the connection 
between the constant parts of the Moon's Parallax and the motions of its Perigee 
and Node. If we had assumed the results of Adams' second theorem (which is 
here a deduction from the equations (A) and (B)), it would have been unneces- 
sary to prove (B). But it is desirable to see exactly by what means the results 
come, and thence to get an idea as to the way in which they might be extended 
to higher powers of e, y. The methods followed above appear to be sufficient to 
indicate further relations if they exist. The chief trouble is to investigate the 
terms depending on both e, y, those depending on one only of these constants 
requiring much less labor. The results are independent of the particular method 
of defining e, y, so long as these quantities are of the same order as the coeffi- 
cients of the principal elliptic term in longitude and of the principal term in lati- 
tude respectively. In fact, their definition is restricted only by the assumption 
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that the coordinates are expressible by means of the same angles at any stage of 
the approximations. We simply approximate to the coefficients. A reference to 
Adams' paper (IV) will make these points clearer. 

17. Extension to Terms of Higher Orders. 

It is not difficult to extend these formulae when terms of higher orders in 
the expressions of c , g are under consideration. But simple relations with the 
constant part of the Parallax of the Moon, like those which exist for q — 1 and 
q = 2, do not seem to be present in the terms of higher orders. The methods 
used here should give any relations which may exist, and it is not at all improb- 
able that there are such. 

The connections between the motions of Perigee and Node and the Parallax 
are really due to the fact that these are all independent (in the sense of the 
theorem of section 11) of the coefficients of the next higher order than the terms 
under consideration in the Perigee and Node. The formulas for the latter, though 
symmetrical in x, y or z, are not so with reference to the quantities x , x e , x y *, 
etc., while the expression for the Parallax is so. Consequently the results of 
Adams' theorems and their extensions must be rather considered from this point 
of view as accidental. 

The terms in the motion of the Perigee depend principally on the function 

(El,_ s +2z 1 z 2q _ z y 

where 2q — 2 is the order of the terms required. 

As an example, take the case g , =3 in the theorem of Art. 11. This gives, 
if we consider only the constant part of the order e 6 , i. e. putting y = 0, 

which may be shown to be equal to 

after the method and notation used before. The first term on the right-hand 
side vanishes, having no coefficient of the order <s 6 ; the second term is that 
required for the determination, in the motion of the Perigee, of the terms of order 
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e 4 [equation (18)], but the other two terms do not vanish. Some reductions are 
possible in these latter terms, but as 4 and y t cannot be eliminated. 

PART III. 

Terms of Long and Short Periods whose Coefficients are Large in 
Comparison to the Orders of the Coefficients. 

There are some results which can be deduced from the infinite determinant 

which are of importance when short-period terms are under discussion. 

Let the arguments of any series of terms in x , y differing by multiples of 

2v (t — t Q ) be 

(2/ -4- A) vt + const., 

where A is of the form j\ + km + pc + %Q9i an( i consequently the coefficients of 
the order a' - 'V*Vy" 9 '. Here k,p,q are integers positive or negative. Also 
j\ z= o when f is even and j\ = 1 when / is odd ; in addition, we have 

&' = (k) or Td = (h) + even positive integer, etc. 

Let the coefficients corresponding to any given value of A be ^, A.J. The 
equations which determine them are 

(2/+ 1 + rn+AfXj +X i M t l J _ i +2,iW-j = known terms independent of %, 7J, 
(2j-l-m+AnL j +2 i M_X-j+2 t y- i ^- i = " " " " 

Suppose that the solution of these equations is found by means of determinants. 
A reference to equations (14) shows that the common denominator in the expres- 
sions of all the coefficients \, Jlj is A (A) or V (A). Now the identity 

V (c) ■= (cos itc — cos nc^cos nc — 1) 

of Art. 4 holds for all values of c, and therefore when we put A for c. Hence, 
dropping the suffix of c x so that c denotes the principal part of the motion of 
Perigee, we have for the common denominator V (A) the value 

(cos 7tA — cos 7tc)(cos 5tA — 1) . 

Prom this we conclude that the coefficients %, %' may become large when either 

cos 7tA — cos tic or cos nh. — 1 
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is a small quantity. Since c is nearly equal to unity, these expressions cannot 
become small together. This shows that when we have the equations 

A = c + even integer or A = even integer, 

approximately satisfied, the quantities %, Jl/ may become large. The term " even 
integer " includes negative as well as positive even integers and also zero. 

The second case is that of long-period terms. The result obtained does not 
add anything to our knowledge of them, except perhaps the fact (which applies 
also to the terms affected by the first equation) that all terms whose arguments 
differ from Avt -\- const, by even multiples of 2v (t — 1 ) should be considered 
together. The first equation affects only short-period terms : the periods of these 
terms differ little from those of the principal elliptic inequalities, that is, from 
2n/(2j± c) v . We have neglected in c all terms but those depending on rn only, 
but the results evidently hold when we consider the full value of c. If we had 
used numerical values for all the constants it would have been necessary to take 
this full value. In the case before us, where we use the numerical value of m 
only, we take the part of c depending on m only. 

It seems to have been usual to consider short-period terms as those whose 
periods approximate to that of the mean motion. The result shows that we 
ought rather to consider the terms whose periods approximate to that of the 
principal elliptic inequality. The result obtained in a previous investigation 
(VII, p. 254) is evidently only a first approximation to that obtained here. The 
difference is due to the fact that here we include at once all the terms whose 
arguments differ from that considered by even multiples of 2v(t — 1 ), and we 
therefore include in the coefficients all powers of m . 

Only a finite number of the coefficients A,-, Xj will become large relatively to 
the rest, but the whole series will be large relatively to their order in e, e 1 , y, 1/a'. 
No coefficients of terms of long period become large by the integration of the 
latitude equation. The same remarks, as to the terms of short period, apply if 
for "principal elliptic inequalities" we substitute "principal latitude inequali- 
ties," i. e. those of period 2it\(2j ± g)v. 

For the calculations of Part I, where we are especially concerned, in the 
infinite determinant, with the values of c, g which depend on m only, the above 
remarks may be included in the following rule : When calculating any series of 
terms in x, y in which A has the value i -f- hm + pc -f- 2qg , in order to discover if 
a large divisor is present it is necessary to examine whether for positive or negative 
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integral values of i, k, p, q either of the expressions 

i-\-km-\-{p — l)c+2qg, i + km -f- pc-\- 2qg 

becomes small. In calculating terms in z, we examine whether the second expression 

may become small. The last statement arises from the fact that in latitude A is 

of the form 

i + km+pc + (2q + l)g, 

and this for large coefficients must be approximately equal to g. It is to be 
remembered that the rule applies to all terms wherein the coefficient of vt in the 
argument differs from A by even multiples of v . 

The part of the coefficient of lowest order of a term inx, y whose argument 
is (i + km -\- pc + 2qg)vt + const., is of order e" e e p y iq or e" c e p y* q .l/a', according 
as * is even or odd: in z the argument being \i-\- Tim-\-pc-\-{2q-\-\))vt-\- const., 
the orders are for i even or odd, e lh e p y 2q+1 and e lk e' p y z ' i+1 l/a' respectively. Hence 
if there is a long-period inequality with a large coefficient of order e" i e p y 2 ' i in x, y 
there is a short-period inequality with a large coefficient of order e" c e F + l y 2q in x, y 
and a similar one of order e" fe e p y 24+1 in z. In the indices of e, e', y the numbers 
k,p, q receive positive values. 

By considering the algebraical values of the principal parts of c, g we can 
also arrive at the order of the denominator in relation to m , and thence, by look- 
ing at the manner in which the denominators are built up, the orders of any 
terms, with respect to the parameter m, can be obtained. 
Oheist's College, Cambridge, December 2ith, 1894. 



